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ABSTRACT 

XTE J1701-407 is a new transient X-ray source discovered on June 8th, 2008. More 
than one month later it showed a rare type of thermonuclear explosion: a long type I 
X-ray burst. We report herein the results of our study of the spectral and flux evolution 
during this burst, as well as the analysis of the outburst in which it took place. We find 
an upper limit on the distance to the source of 6.1 kpc by considering the maximum 
luminosity reached by the burst. We measure a total fluence of 3.5xl0 -6 erg/cm 2 
throughout the ~20 minutes burst duration and a fluence of 2.6 xlO -3 erg/cm 2 during 
the first two months of the outburst. We show that the flux decay is best fitted by a 
power law (index ~1.6) along the tail of the burst. Finally, we discuss the implications 
of the long burst properties, and the presence of a second and shorter burst detected 
by Swift ten days later, for the composition of the accreted material and the heating 
of the burning layer. 
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1 INTRODUCTION 

The bulge of our Galaxy harbours a large fraction of the 
known population of accreting compact objects. During 

' one of the regular monitoring observations of that region 
with the proportional counter array (PCA) onboard the 
Rossi X-ray timing explorer (RXTE) on June 8th, 2008, 
a ne w transient X-ray source was discovered: XTE J1701- 

' 407 l|Markwardt et all l2008ah . Three days later a follow- 
up observation with the X-ray telescope (XRT) onboard 
Swift improved the source localization and showed a spec- 
trum that could be fitted with an abs orbed power law of 
index ~2 (De genaar fc Wiinandsl [2008h The transient na- 
ture, flux level and spectrum of the source are character- 
istic of X-ray binaries, but the exact class of the binary 
(high-mass or low-mass) was unkown and the compact ob- 
ject (black hole or neutron star) remained unidentified for 
more than one month. On July 17th at 13:29:59 UT the 
burst alert telescope (BAT) onboard Swift detected an X- 
ray flare at a position co nsistent with that of XTE J1701-407 
(Bart helmv et al1l2008l h Swift began slewing to the source 
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60 seconds after the trigger and the XRT began observing 
the field 97 seconds after the initial BAT trigger, detect- 
ing a fading X-ray source a t the position of XTE J 1701- 
407 (|Barthelmv et alj|2008h . iMarkwardt et al.l (|2008bh sug- 
gested that the flare could be caused by a thermonuclear 
burst based on the early BAT data. By study ing the evolu- 
tion of the XRT spectrum. lLinares et al . ( 2008) found a clear 
cooling curve and confirmed the speculation that this was 
a thermonuclear event, identifying the source as an accret- 
ing neutron star (NS) and the system, in all li kelihood, as a 
low-m ass X-ray binary (LMXB). As noted bv lLinares et al.1 
(2008) this particular event, having a total duration of more 
than 15 mi nutes, belongs to the r are subclass of long dura- 
tion bursts (|Cumming et al.ll2006l ). 

After the long X-ray burst the source continued 
in outburst and R XTE kept observ i ng it . A refined 
Swift-XRT position l|Starling fc Evansl l2008h . a search 



1 To clarify the terminology we stress that in the context of 
LMXBs an outburst is powered by accretion and lasts for weeks 
to years whereas a type I X-ray burst, or "burst", results from 
unstable thermonuclear burning on the surface of a NS and lasts 
for ~10 seconds to ~hours. 



© 0000 RAS 



2 Linares et al. 



for an IR counterpart |Kaplan fc Chakrabartvl I2008T I and 
the discovery of kHz quasi-periodic oscillations (QPOs; 
IStrohmaver et all 120081 ) were all reported during the 
following two weeks. On July 27th BA T detected an- 
other X-ray flare from XTE J 170 1-407 ^Sakamoto et all 
2008). The soft spectrum and duration (~10 s) of 
the flare indicated that this was a type I X-ray burst 
| |http:// gcn.gsfc.nasa.gov/notices-s/318166/BA/[), and 
hereafter we refer to it as the short burst (even though from 
the BAT data the total duration is somewhat uncertain). 
At the moment of writing the source is still active, more 
than two months after the start of the outburst. 

In this Letter we present our analysis of the bursting 
behaviour of this source. Long bursts are rare, and there- 
fore probe quite unusual burning regimes. Their duration 
is between that of normal Type I X-ray bursts (durations 
~ 10 — 100 s, triggered by unstable burning of H/He) and su- 
perbursts (durations ~ hours, triggered by unstable burning 
of C) . There are two scenarios in which such bursts could oc- 
cur, both requiring the build-up and ignition of a thick layer 
of He. One possibility is that th e system is ultracompa ct, so 
that it accretes nearly pure He l|Cumming et al. 2006). The 
other possibility is th at unstable H bur n ing at low accretion 
rates builds up He (|Peng et al.l 120071 ; ICooper fc Naravanl 
2007). While several e xamples of the first type are now 
know n (see for example lin't Zand et alj|2007l ; iFalanga et al.l 
120081 and references therein), the second class - which probe 
the boundar ies of H burning stabi lity - have so far proved 
much rarer |Chenevez et al.l 12007), In this Letter we will 
demonstrate that XTE J1701-407 is a good candidate for 
membership of this second class. 



at 13:31:36.9 UT (Sec. [TJ and lasted for a total of about 
one kilosecond. All the XRT data were collected in win- 
dowed timing (WT) mode, except for the initial ~8 seconds 
which were taken in photon counting (PC) mode. We di- 
vided the observation into nine contiguous intervals of in- 
creasing length in order to compensate for the flux decay 
and maintain an approximately constant total number of 
counts in each interval (between ~1500 and ~2500). About 
830 seconds after the start of the observation a data gap of 
~3300 seconds occurred, after which the XRT continued to 
observe the source for another ~35 seconds. This last data 
segment allows us to estimate the persistent flux after the 
X-ray burst. We extracted source and background spectra 
for each interval using circular regions with radii of ~30 and 
~15 pixels, respectively, and fitted the resulting spectra us- 
ing an absorbed black body model (reduced \ 2 between 0.7 
and 1.1). 

We performed fast Fourier transforms on 2s segments 
of XRT-WT data taken during the long burst in order to 
search for rapid variability. For that purpose we used the 
energy range 0.5-10 keV, keeping the original time resolution 
and thereby sampling the ~0. 5-280 Hz frequency range. 

Both post-burst observations (00317205001 and 
00318166000) allow us to characterize the persistent soft 
X-ray emission of XTE J1701-407. For the PC data we 
used an annular extraction region with radii 10-25 pixels in 
order to account for pile-up. The energy fraction enclosed 
by the resulting annuli was about 20%, and such energy loss 
was taken into account when creating the ancilliary files. 
Finally, including also the WT data, we fitted the energy 
spectra with an absorbed power law model (reduced x 2 of 
1.1 and 1.2). 



2 OBSERVATIONS AND DATA ANALYSIS 

We analysed all the 27 pointed RXTE-PCA observations of 
XTE J1701-407 taken until August 10th, 2008 (ID 93444- 
01). We also analysed the Swift-'KKT data taken during 
and after the long burst (on July 17th, obsids 00317205000 
and 00317205001). Swift observed XTE J1701-407 again 
on July 27th, triggered by the BAT detection of the short 
burst (Sec.[T]). We also included this XRT observation (obsid 
00318166000) in our analysis, which started ~109 seconds 
after the trigger but did not detect the burst. In the follow- 
ing two sections we provide the details of our analysis. All 
errors given in Section correspond to the 1 a confidence 
level. 



2.1 Swift 

We ran the XRT pipeline (v. 0.12.0) on all the Swift obser- 
vations. For the spectral analysis we created exposure maps 
and ancilliary files and used the latest response matrices 
(vOll) provided by the Swift-XBT team. We grouped the 
spectra to a minimum of 20 counts per bin, applyied a sys- 
tematic error of 2.5% l|Campana et al. 20081) and fitted the 
0.5-10 keV spectra within Xspec (|Arnaudlll996l . v. 11). In 
our fits (see models below) we fixed the absorbing column 
density to the average value found from X-ray fits to the 
persi s tent emission: 3.5xl0 22 cm -2 |Degenaar fc Wiinands! 
120081 : iMarkwardt et al.ll2008al ). 

Observation 00317205000 started on July 17th, 2008 



2.2 RXTE 

We analysed all the pointed RXTE observations of the 
source made until August 10th. After applying standard fil- 
ters we extracted energy spectra from Standard 2 data col- 
lected by PCU2 and corrected them for background using 
the latest PCA background models. We generated reponsc 
matrices and fitted the spectra within Xspec in the 3.0-25.0 
keV band, after applying a 1% systematic error, with a black 
body plus power law model corrected for absorption (see 
Sec. 12. ip . The resulting fits had reduced x 2 between 0.5 and 
1.4 (for 45 degrees of freedom). We extracted a hardness ra- 
tio ("hard colour" hereafter) for each observation, defined 
as the count rate in the 9.7-16.0 keV band divided by that 
in the 6.0-9.7 keV band and normalized by that of the Crab. 
Furthermore, we extracted a background and deadtime cor- 
rected HEXTE (cluster B) energy spectrum from each ob- 
servation, but found no significant source emission above 50 
keV in any of them. 



3 RESULTS 
3.1 Outburst 

The overall lightcurve is shown in Figure[T] together with the 
evolution of the hard colour. The outburst shows two peaks 
at a similar flux, (6.2, 6.8) x 10" 10 erg/cm 2 /s. The BAT tran- 
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Figure 1. Time evolution of the 2-50 keV unabsorbed flux (top, 
filled dots) and the hard colour (bottom) during the outburst of 
XTE J1701-407, as measured by the RXTE-PCA. Arrows indi- 
cate the times of the long (LB) and short (B) X-ray burst. The 
15-50 keV flux measured by Swift-HAT is shown as well (top, grey 
rectangles). 

sient monitor lightcurveijf] indicate that the source was ac- 
tive during the three weeks where no RXTE pointed observa- 
tions were obtained, at a flux level similar to that measured 
during the first two weeks of outburst (see Fig. [TJ . The hard 
colour traces directly the changes in spectral state along the 
outburst, with a clear tendency: when the source is bright 
its spectrum is soft whereas the spectrum is harder at the 
lowest fluxes. However, it is interesting to note that as is the 
case for many other NS-LMXBs the correspondence bet ween 
flux and hardness is not one-to-one (|van der Klij |2006V To- 
gether with the presenc e of kHz quasi-periodic oscillations 
l)Strohmaver et aTll2008t ). the peak luminosity and the spec- 
tral evolution (Fig. [TJ show that XTE J1701-407 is a new 
mem ber of the atoll source class (jHasinger fc van' der Klisl 
1989). Both X-ray bursts occur at relatively high fluxes (see 
Table [TJ and soft spectra (low hard colour). 

By interpolating linearly between the available data- 
points we integrate the lightcurve and find a total outburst 
fluence up until August 10th of 2.6xlCP 3 erg/cm 2 . More 
than half of this energy (1.6xl0 -3 erg/cm 2 ) was radiated 
before the long burst occurred, whereas the fluence between 
long and short bursts was about 5xl0 -4 erg/cm 2 . 

From our spectral fits to the post-burst XRT obser- 
vations (Sec. 12. ip we measure 2-10 keV fluxes of (4.0 ± 
0.9) x 10~ 10 erg/cm 2 /s (July 17th) and (5.6 ± TO) x 10~ 10 
erg/cm 2 /s (July 27th), which are consistent with the fluxes 
in the same band obtained from the nearest (within one 
day) PCA observations. In both cases the photon index was 
~2. Furthermore, we estimate unabsorbed 1-2 keV fluxes of 
2.3 x 10" 10 erg/cm 2 /s (July 17th) and 3.5 x 10~ 10 crg/cm 2 /s 
(July 27th) and add them to the unabsorbed 2-50 keV 
fluxes measured with the PCA within one day of the bursts 
(6.3 x 10" 10 erg/cm 2 /s for both bursts ; see Fig. [1} . These 
1-50 keV fluxes (Table[TJ are the closest we can safely get to 
the bolometric persistent flux at the times of the bursts, but 

2 http:/ /swift. gsfc.nasa.gov/docs/swift/results/transients/ 



we note that the uncertainties in the absorbing column den- 
sity and the extrapolation to the 0.01-1 keV energy range 
could change the bolometric flux b y a factor of ~2 (see e.g. 
discussion in [in 't Zand et al]|2007h . 



3.2 X-ray bursts 

From the flux and sp ectral parameters reported by 
iMarkwardt et all (|2008bl ) we estimate a bolometric unab- 
sorbed peak flux of 8.4xl0 _s erg/cm 2 /s for the long burst. 
By assuming that the peak luminosity was equal to or 
lower than that typical of phot ospheric expans i on ty pe I 
X-ray bursts (~3.8xl0 38 erg/s iKuulkers et alJ [20031 ) we 
can place an upper limit on the distance to the source 
of 6.1 kpc. On the other hand, if we use the Edding- 
ton limit for a mix of H and He (~1.6xl0 38 erg/s for a 
hydrogen fraction X=0.7) as a limit for the peak lumi- 
nosity the upper limit on the distance becomes 4.0 kpc. 
The rise time in the BAT lightcurve was about 50 sec- 
onds, the peak occurred about ten seconds after the trigger 
(http://gcn.gsfc.nasa.gov/notices_s/317205/BA/) and the 
total duration, assuming that the burst finished when it 
reached 1% of the peak flux, was about 21 min. The fluxes 
in the post-burst XRT observation, taken more than two 
hours after the peak of the long burst, are similar to the 
flux we measure in the last interval (~1 hour after the peak, 
Sec. 12. 1[) , confirming that the source had returned to the 
persistent level. 

Figure [2] shows the results of our time-resolved spec- 
troscopy of the tail of the long burst. The flux decays by 
almost two orders of magnitude and the temperature de- 
creases from 1.72± 0.05 keV to 0.88±0 .02 keV during the 
burst tail (see also iLinares et al1l2008h . XRT did not ob- 
serve the peak of the burst, and n o spectral variations wer e 
seen in BAT data of that period (|Markwardt et al.ll2~008bl ). 
Therefore the question of whether or not the burst showed 
photospheric radius expansion cannot be addressed. Going 
back to our XRT measurements, the radius of the black 
body appears to increase between ~200 s and ~500 s af- 
ter the peak (Fig. [2}. However, the large error affecting this 
parameter, the uncertainties on the underlying model (NS 
atmosphere in a dramatically-non-steady state) and the ab- 
sence of an associated drop in the black body temperature 
imply that this cannot be taken as concrete evidence of a 
physical increase in emitting area. 

We fit the flux decay of the long burst measured by 
the XRT and find that a simple power law with an index 
of -1.62T0.04 and normalization (3.0T0.6) x 10~ 5 erg/cm 2 /s 
gives a remarkably good description of the data (reduced 
X 2 of 1.2; see Fig. [3|. A simple or double exponential de- 
cay do not fit the data satisfactorily (reduced \ 2 °f 22 and 
7, respectively). We note that power law cooling is con- 
sistent with the expectation s from models for long bursts 
|Cumming fc Macbeth! |2004| ). 

The bolometric fluence measured by XRT during the 
long burst, after subtracting the persistent flux level, is 
2.2 xl0~ 6 erg/cm 2 . We integrate the flux decay along the 
data gap (Sec. 12. ip and estimate a fluence of 4.4xl0~ 7 
erg/c m 2 during that period. A dding the reported BAT flu- 
ence l|Markwardt et alj|2008bh we find a total burst fluence 
of 3.5xl0 -6 erg/cm 2 , which using the 6.1 kpc upper limit 
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Figure 2. Time evolution of the bolometric flux (top), black 
body temperature (middle) and black body radius (bottom) along 
the long X-ray burst. Filled dots correspond to our Swift-XIIT 
measurements whereas e mpty squares show Swift-BAT results 
llMarkwardt et al"1l2008bT ), The horizontal dotted line on the top 
panel indicates the persistent flux. Error bars in the x axis repre- 
sent the averaged intervals. The distance upper limit of 6.1 kpc 
has been used to obtain the black body radius. 



on the distance corresponds to a maximum radiated energy 
of 1.6xl0 40 erg. 

We searched the XRT power spectra obtained during 
the tail of the long burst (Sec. 1 2 . 1 p and found no strong 
pulsations or QPOs. This non-detection is not surprising 
given i) the distrib ution of burst oscillation frequencies 
l|Gallowav et alj2008ft combined with our Nyquist frequency 
(~280 Hz) and ii) the relatively low count rate collected by 
the XRT: for the maximum rate of ~90 c/s, a 3<r detection 
of a 2 Hz wide QPO i n a 2s FFT would require a fractional 
rms of at least ~ 25% (Ivan der Klidl 1995ft, much higher than 
the values usually measured ( Galloway et alj|2008ft . 



4 DISCUSSION 

Theoretical modelling of long bursts suggests that they in- 
volve the build-up and ignition of a thick layer of He. One 
possibility is that the binary is ultracompact, so that the 
accreted fuel is almost pure He. At low accretion rates, 
and in the absence of heating from steady H burning, the 
temperature remains low and ignitio n can be delayed un- 
til a thick layer of H e has built up (lin't Zand et al.ll2005l ; 
ICumming et alj|2006ft . Although He burning is quick, the 
cooling time is long because of the thickness of the layer 
in which the heat is deposited. Tong bursts can also be 
triggered in systems accreting a mix of H and He if the 
source accretes at low rates, close to the point where H 
burning stabilises ( Fuiimoto et alj Il98ll ; iPeng et af 1 120071 : 
ICooper fc Naravanl 2007). In this case weak H flashes can 
build up He, which ignites sporadically, resulting in a long 
burst. Layer thickness again results in a long cooling time: 
but the presence of H in the burning mix may also prolong 
nuclear energy generation. 
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Figure 3. Fit to the long burst flux decay with a power law of 
index ~1.6, together with the values of the baseline-subtracted 
bolometric flux. 



Using the measured fluxes F from Table [T] and the up- 
per limit on the distance d that we obtained from the long 
burst, we infer an upper limit on the global accretion rate 
M w 4nd 2 RF/GM of 3.7xl0~ 10 M Q yr" 1 (for a neutron 
star mass M = IAMq and radius R — 10 km) around 
the times of the two bursts. Assuming isotropy, this cor- 
responds to an upper limit on the local accretion rate m of 
1.9xl0 3 g/cm 2 /s, (2.5% of the Eddington rate for an Ed- 
dington luminosity of 1.6xl0 38 erg/s). This accretion rate 
is close to the rate where H burn i ng is expected to tran- 
sition to stability jFuiimoto et al] Il98ll : IPeng et al.l 120071 ; 
ICooper fc Naravanl l2007), so we need to examine both sce- 
narios for the generation of the long burst. 

The BAT rise time of the long burst, ~ 50 s, is 
mu ch longer than most of th e pure He bur s ts exa mined 
by lin't Zand etafl (120071 ) and iFalanga et all [|200Sft . This 
points to it being a mixed H/He type burst, triggered by H 
ignition, although we note that burst lightcurves can change 
considerably depending on the energy range where they are 
observed (see IChelovekov et al.ll200fj : lMolkov et alj|2005ft . 

We next consider the ignition conditions. The long burst 
had a total energy output of at most Et = 1.6 x 10 40 ergs 
(6.1 kpc distance). The corresponding ignition column depth 
y — Eh(l + 2)/47ri? 2 Q nuc , where z is the redshift, R the neu- 
tron star radius and Q lluc = 1.6 + 4X MeV/nucleon the 
nuclear energy releas e, given an H fraction X at ignition 
(|Gallowav et al.ll200Sft . Assuming z = 0.31 (the redshift for 
a 1.4 Mq, 10 km radius neutron star) we derive an ignition 
depth y = 3.9 x 10 s g/cm 2 for solar abundances (X=0.7) and 
y — 1.1 x 10 9 g/cm 2 for pure He (X=0). Ignition of pure He 
at this column depth would require a heat flux from the 
deep crust of > 2 MeV/nucleon at 1% Eddington accretion 
rates (Figure 22, ICumming et al.ll2006ft . This is higher than 
th e crustal energy rel e ase of ~ 1.5 MeV/nucleon calculated 
bv lHaensel fc Zdunikl <|l990l , l2003ft . but close to the value of 
~ 1.9 MeV/nucleon rece ntly found by iGupta et ail [|2007ft; 
lHaensel fc Zdunikl | 2008ft . Recent results bv lHorowitz et all 
( 2008ft demonstrate that even larger amounts of heat may 
be released in the crust of an accreting neutron star. How- 
ever, the high heat requirements certainly put some strong 
constraints on the pure He accretion scenario. Using the 
maximum local accretion rate calculated above (1.9xl0 3 
g/cm 2 /s), these ignition conditions predict minimum recur- 
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Table 1. Properties of both type I X-ray bursts detected from 
XTE J1701-407 



Lone 



Short 



BAT rise time (s) 

Duration (min) 

Peak flux^ (10 -8 erg/s/cm 2 ) 

Total fiucncc^ (10~ 6 erg/cm 2 ) 

Persistent flux£ (10~ 10 erg/s/cm 2 ) 



50 
21 
8.4 
3.5 
8.6±0.9 



5 
? 

~4.9 

? 

9.8±1.0 



a From start of rise to 1% of peak flux. 

^ Bolometric 

c Including fluence in BAT l|Markwardt et al.ll2008bh and 
during XRT data gap. 
d Unabsorbed 1-50 keV 



rence times At If!C = y(l + z) /rh of ~10 days (X = 0) and ~3 
days (X = 0.7). The data provide only weak constraints on 
the time without bursts that elapsed before the long burst 
(we estimate that the source was in the field of view of BAT 
for about 20% of the two months analysed herein), but these 
numbers are at least compatible with the time since the start 
of the outburst. 

The presence of a short, weaker burst at similar accre- 
tion rates is perhaps the strongest piece of evidence point- 
ing to mixed H/He fuel. The short burst has a relatively 
slow rise time, « 5 s, similar to that expected for mixed 
H/He bursts. The fact that it has a lower peak flux than 
the long burst also implies some H content: if the system 
were a pure He accretor we would expect both bursts to 
reach similar peak fluxes since He bursts, whether long or 
short, are expect ed to exhibit photospheric radius expansion 
ICummindbOoj ). In the H-triggered burst scenario it is pos- 
sible to have bursts with quite different properties at very 
similar accretion rates, as H burning transitions to stability 
l|Cooper fe Naravanll2007l ). 

XTE J1701-407 increases the number of known NS- 
LMXBs in the interesting low accretion rate bursting regime. 
The balance of evidence suggests that it may be a good 
probe of thermonuclear burning near the boundary of sta- 
ble H ignition. However, further observations, in particular 
better constraints on recurrence times and burst energet- 
ics, are required to confirm this picture. Identification of the 
companion star and determination of the orbital parameters 
may also be able to confirm whether or not the system is 
ultracompact. 
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